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Abstract 

The ALICE experiment has measured quarkonia production in pp and Pb-Pb collisions at the 
CERN LHC, in the rapidity ranges |y| < 0.9 and 2.5 < y < 4. Quarkonia are considered to 
be a sensitive probe of deconfinement, and a detailed differential study of their yields can give 
important information on the properties of the medium created in heavy-ion collisions. In this 
paper, we will mainly discuss the centrality dependence of the J /if/ nuclear modification factors, as 
well as their pj and y dependence in bins of centrality, which will be then compared to theoretical 
models. Preliminary results on the J /if/ elliptic flow and on if/(2S ) production will also be shown. 



1. Introduction 

The suppression of quarkonium states in heavy-ion collisions was proposed long ago JJ] 
as a clean signature of a formation of a Quark-Gluon Plasma (QGP). In addition to the color 
screening mechanism proposed in [1] as the main source of suppression, it became soon clear 
that other effects may contribute, as dissociation in cold nuclear matter |2[] and/or in the hot 
confined medium [3]. Also initial state effects, such as nuclear shadowing and initial state parton 
energy loss are expected to play a role |0]. On the experimental side, lots of data are now 
available at SPS [5, 6] and RHIC energies, and a quarkonium (Q) suppression beyond cold 
nuclear matter effects (the so-called "anomalous suppression") was detected, either by studying 
the ratio ctq j crDreii-Yan or the nuclear modification factor Raa, for different collision systems and 
as a function of the centrality and of the quarkonium kinematic variables |9J]. It is commonly 
accepted that the observed anomalous suppression is likely to be related to the formation of a 
deconfined state in nuclear collisions. In this situation, the first ALICE results [ 10] brought in a 
considerable surprise, showing for the J /if/ clear hints of a smaller suppression, at low transverse 



momentum, compared to observations at lower energies. Partonic transport models 11UI12I1 . as 
well as statistical generation models 11311 . interpret such an effect as being due to a (re)generation 
of J/if/ along the collision history and/or at hadronization, favoured by the large cc multiplicity 
occurring at LHC energy. A confirmation of this interpretation clearly requires a comprehensive 
study of the i/i// yield in nuclear collisions, with an emphasis on its dependence on kinematic 
variables. Such a study requires a sizeable Pb-Pb statistics, which was indeed collected by the 
ALICE experiment at the end of 2011, and results from this data taking are the main object of 
this paper. An accurate set of pp data is a fundamental pre-requisite for a study of nuclear effects 
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on any observable. Moreover, important QCD-related topics can be investigated by the study of 
quarkonium hadro-production in elementary collisions. We will therefore also summarize the 
main outcomes of ALICE pp studies on quarkonia. 

2. ALICE experiment and data taking conditions 

ALICE is the dedicated heavy-ion experiment at the CERN LHC. Details on the experimental 
set-up can be found in [14]. The measurement of quarkonia production is carried out in the 
central barrel (|y| < 0.9) through their e + e~ decay, while at forward rapidity (2.5 < y < 4$ the 
ptpT decay is studied in the muon spectrometer. 

For the measurement in the central barrel, the data acquisition is triggered by a minimum bias 
condition (MB), essentially defined using information from the two innermost planes of the Inner 
Tracker System and from two scintillator arrays (VZERO) placed at forward rapidities (-3.7< r\ 
<-1.7, 2.8< J] <5.1). In Pb-Pb, several centrality classes can be further defined at the trigger 
level by means of thresholds on the total signal amplitude in the VZERO. For the measurement 
in the forward muon spectrometer, one can require, in addition to the MB condition, one or two 
(depending on the data taking) candidate muons to be detected by the muon triggering system. 

The analysis presented in this paper is based, for Pb-Pb collisions at central rapidity, on a 
combination of 2010 data, taken with a MB trigger, and 201 1 data where, in addition, centrality 
selections corresponding to very central events (~ 0- 10%) and semi-central events (~ 10-50%) 
were performed. The integrated luminosity is L mt ~ 15 fib~ l . The Pb-Pb forward rapidity results 
refer to 201 1 data (L; nt ~ 70 fib ), and were collected by requiring a dimuon trigger, which can 
select both opposite- and like-sign pairs. In order to properly define Raa, one needs a sample 
of pp data collected at the same energy (^Js = 2.76 TeV) and in the same kinematic domain of 
Pb-Pb data. Such data were taken in 2011, collecting L; nt ~ 1.1 nb~' with a MB trigger for the 
central rapidity analysis and L; nt ~ 20nb~' with a single-muon trigger for the forward rapidity 
analysis. Finally, pp data at -\fs = 7 TeV, the top energy of the 2010/201 1 LHC data taking, were 
also collected and analysed, the results presented here referring to integrated luminosities up to 
-100 nb 1 . 

3. Data analysis 

The signal extraction is performed by means of an analysis of the invariant mass spectra, for 
both e + e~ and fi + fT decay modes. In the electron channel a large combinatorial background is 
present, and its size is evaluated using an event mixing technique. The mixed event sample is then 
normalized to the opposite sign (OS) spectrum in the mass region beyond the J/i/f and subtracted 
out. The number of signal events is obtained by bin counting in the mass range 2.92-3. 16 GeV/c 2 , 
the factor extrapolating to the full mass range being obtained via Monte Carlo (MC) simulation. 
Fig.Q](left) shows the invariant mass spectra before and after background subtraction for the case 
where the latter is larger, i.e. central Pb-Pb collisions. 

In the muon channel, the J/if/ signal is extracted by means of a fit to the OS mass spectrum, 
using a phenomenological shape for the background, and a Crystal Ball shape for the signal. The 
result, for centrality integrated Pb-Pb collisions is shown in Fig.[T](right), and the total statistics 



2 In the ALICE reference frame, the muon spectrometer covers a negative pseudorapidity interval (-4 < J] < -2.5). 
Since in pp and Pb-Pb the physics is symmetric with respect to y=0, we drop the negative sign when quoting rapidity. 
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Figure 1 : Left: invariant mass spectrum for Pb-Pb (0-10% central) opposite-sign electron pairs and background calculated 
using a mixed event technique (top), background subtracted spectrum with MC expected shape superimposed (bottom). 
Right: invariant mass spectrum. The result of the fit around the J/i/i invariant mass region is superimposed. 

amount to ~ 4 • 10 4 signal events (the if/(2S) signal contribution has no practical effects on the 
J/if/ yield, and was neglected in this fit). 

The product A x e of acceptance times efficiency was calculated via MC using as input sim- 
ulated pp or Pb-Pb events. Typical values amount to ~ 5 - 10% in the \y\ < 0.9 region (elec- 
trons) B15I1 and 15 - 30% for 2.5 < y < 4 (muons), the latter range of values mainly depending on 
the choice of the trigger threshold and on the request for 1 or 2 matched muons between tracking 
and trigger. For the muon channel in Pb-Pb, the calculation was done embedding MC J/i// inside 
real MB collisions. The centrality dependence of A x e in Pb-Pb collisions is found to be weak, 
not exceeding ~ 8% between central and peripheral events in the muon analysis, and even less 
in the electron analysis. 

For pp collisions, the J /if/ yields are then normalized to the integrated luminosity, calculated 
starting from the number of collected MB events. The absolute J /if/ cross sections are then ob- 
tained relatively to the MB cross section, obtained via a van der Meer scan. The total systematic 
uncertainty on the cross section measurements ranges, depending on the data taking, between 8 
and 12% in the muon channel and between 14 and 18% in the electron channel. 

For Pb-Pb collisions, the results are given in terms of the nuclear modification factor Raa, 
calculated normalizing the J /if/ yield to the corresponding pp cross section times the nuclear 
overlap function, derived from a Glauber calculation. On these analyses, an important source of 
(global) systematic uncertainty is related to the total error on the pp reference cross section. It is 
by far dominant for the electron analysis (~26%) and represents a sizeable contribution also in 
the muon analysis where for example it accounts for about 50% of the total uncertainty for the 
integrated Raa- 

4. Results 

We first briefly summarize the main results obtained by ALICE in the study of pp collisions. 
These include the calculation of the inclusive production cross sections at -\fs — 7 and 2.76 TeV. 
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At 7 TeV the differential cross sections d 2 cr/dydpj were obtained at both central and forward 
rapi dity 111 611 . where they showed a good agreement with results from the other LHC experiments 
(see IU6I1 for a comparison) in the common kinematic ranges, and were unique for what concerns 
the low-px region at midrapidity. At -\fs — 2.76 TeV a differential study was possible only in the 
forward region, while at midrapidity a /^-integrated cross section was obtained, due to the rather 
small MB integrated luminosity II 1711 . 

Other interesting studies performed in pp include the first study at LHC energy of the in- 
clusive J /if/ polarization. The results, published in [18], show that J /if/ production is essentially 
unpolarized up to pj — 8 Ge V/c. The study of this observable has critical implications for the the- 
oretical understanding of quarkonium production, since models were quite unable to reproduce 
coherently the results obtained at lower energy. However, first comparisons of ALICE results 
with NLO NRQCD calculations show an encouraging agreement in the explored pj range fl9ll . 
Another study unique to ALICE is the determination of the dependence of the inclusive J/if/ yield 
on the charged multiplicity measured at central rapidity. The results, published in [20], show a 
linear increase of the yield with charged multiplicity, an observation still waiting for a satisfac- 
tory explanation from theory. Finally, at midrapidity, ALICE has provided a first measurement 
of the fraction of ]/if/ from b-decays at low pj, which also leads to an estimate of the total bb 
cross section 02111 . 



1.2 

1 

0.8 
0.6 
0.4 
0.2 




ALICE, Pb-Pb ^ = 2.76 TeV 
|y|<0.9, p >0 GeV/c, L.-15 jib 1 

PHENIX (PRL 98 (2007) 232301 ), Au-Au^ = 0.2 TeV 
|y|<0.35, p >0 GeV/c global sys.= ± 1 2% 



100 150 200 250 300 350 400 

(N > 
x part' 



1.2 



0.6 
0.4 
0.2 



ALICE Pb-Pb ^ = 2.76 TeV 
■ 2.5<y<4, p >0, L = 70 [lb' 1 

PHENIX (PRC 84(2011) 054912), Au-Au ^ = 0.2 TeV 

□ Inclusive J/i|/, 1.2<|y|<2.2,p >0 GeV/c global sys.= ± 9.2°/ % 



50 100 



150 200 250 300 350 400 
(N > 

part 



Figure 2: Inclusive J/1/7 Raa vs the number of participant nucleons N pan , compared with results from PHENIX. The left 
plot corresponds to the central rapidity results, the right plot to forward rapidity. The boxes around the points represent 
the uncorrelated systematic uncertainties, while the filled boxes on the right correspond to global systematic uncertainties. 



Moving to Pb-Pb collisions, we present in Fig.|2]the centrality dependence of the inclusive J/if/ 
nuclear modification factor (see also 11221,1231, 1 1510 . compared with results obtained by PHENIX. 
While lower energy results show an increasing suppression moving towards central collisions, 
ALICE results clearly indicate a saturation of the suppression at both central and forward rapidity. 
While no final conclusion can be drawn on the relative size of the suppression measured in 
ALICE in the two rapidity ranges, mainly due to the large normalization error corresponding 
to the reference pp cross section, there is a clear evidence for a smaller suppression at LHC 
with respect to RHIC energy. Partonic transport models which include a (re)generation process 
for J /if/ due to the recombination of cc pairs along the history of the collisions indeed predict 
such a behaviour ||3lllllll2|], the smaller suppression at the LHC being due to the larger cc pair 
multiplicity which compensates the suppression from color screening in the deconfined phase. 
A similar behaviour is obtained by the statistical model lfl3ll . where the J /if/ yield is completely 
determined by the chemical freeze-out conditions and by the abundance of cc pairs. 
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A deeper insight of the mechanisms at play in Pb-Pb at the LHC can be obtained by studying 
the centrality dependence of Raa in bins of transverse momentum. If indeed (re)generation 
mechanisms play a role, their impact should be much more important, for simple kinematic 
arguments, at low transverse momentum. A larger Raa should then be observed, in particular for 
central collisions, at low pj with respect to high pj. In Fig. [3] we separately show the centrality 
dependence of the J /iff Raa for < pi < 2 and 5 < pj < 8 GeV/c. A larger Raa can indeed 
be observed at lower pj. The results are compared with the calculation of a partonic transport 
model [11] which is in good agreement with the data, in particular for central collisions. The 
contribution of (re)generated J/tfr is shown as continuous lines in the plots and indeed is sizeable 
only at low p^. Similar conclusions can be reached bya direct study of the pi dependence of 
Raa-, both centrality integrated and in bins of iV part B22l 12311 . 
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Figure 3: Raa vs centrality for two transverse momentum bins (0< pj <2 and 5< pj <8 GeV/c), superimposed with the 
results of a partonic transport model. The contribution of primordial (dashed) and regeneration (continuous line) compo- 
nents is shown, as well as their sum. The area between the curves represent an estimate of the theoretical uncertainty. 



Other interesting information can be obtained by studying the rapidity dependence of Raa- 
RHIC results showed a significantly larger suppression at forward rapidity [7], an observa- 
tion which lead to many theoretical speculations, mainly connected with the observation of a 
(re)combination component at central rapidity, where the charmed pair density is larger, or with 
a different size of cold nuclear matter effects |2D. In Fig.|4](left) we show the centrality integrated 
Raa vs rapidity. We note that the central rapidity result comes from the lower statistics 2010 Pb- 
Pb run, since in the high luminosity 2011 data taking the bulk of the statistics was mainly for 
central events, due to the trigger selection. One can see, within the forward rapidity range, an 
increase of the suppression moving to larger y. In Fig. [4] (right) we present the centrality de- 
pendence of Raa for various rapidity ranges. The tendence to a larger suppression at forward y 
can be seen at all centralities. Again, the large normalization error (filled boxes) prevents a firm 
conclusion from the comparison between central (electron) and forward (muon) results. It can 
be noted that shadowing calculations performed using recent parameterizations (EKS09, nDSg) 
suggest that the suppression observed at midrapidity could be compatible with a pure nuclear 
shadowing scenario, while other effects must be invoked in order to reproduce the Raa observed 
at forward rapidity EH]. 

Finally, it must be noted that for the moment no data on cold nuclear matter effects on J/i// 
exist, at LHC energy. A precise determination of shadowing effects and of the J/i// break-up 
probability via interactions with cold nuclear matter (the latter effect might become negligible at 
the LHC, due to crossing times much smaller than charmonium formation time) can be obtained 
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Figure 4: Left: i?AA vs y for Pb-Pb collisions. Global systematic uncertainties are shown in the legends. Right: Raa 
vs centrality in y bins. For the central rapidity points, the boxes show the uncorrelated systematic uncertainties. For the 
forward rapidity points the shaded boxes represent uncorrelated systematic uncertainties, while the brackets correspond 
to uncertainties correlated for that rapidity bin, but uncorrelated between various rapidity bins. Finally, the filled boxes 
on the right represent global systematic uncertainties, uncorrelated between central and forward rapidity. 



by studying p-Pb collisions. Such a run will take place at the beginning of 2013 and will help to 
make the current conclusions from Raa studies more robust. 

The differential Raa results previously discussed can be considered as strong hints for an 
important contribution of (re)generated ]/iff at low pj in Pb-Pb collisions at the LHC. An in- 
dependent confirmation of this hypothesis may come from the study of the elliptic flow of J/ifr. 
If (re)generation effects are sizeable, then the corresponding i/tfr would inherit the flow related 
to the collective expansion, which is experienced by the charm quarks contained in the fireball. 
Theoretical models predict in this case a non-zero V2 for the ]/if/ at intermediate pi ll24ll . 

For this study, which was performed on the forward rapidity event sample, methods based 
on the event plane determination have been used. The signal in the VZERO detector was used 
for such a purpose, and the usual fit of the J/if/ azimuthal distribution adopted for these studies, 
diVW/dA^ = A X (1 + 2v° bs cos 2A0) was performed, where A0 is the azimuthal angle of the J/ift 
measured with respect to the event plane. The v^ s quantity is then corrected for the event plane 
resolution, and plotted in Fig.|5](left) in four pj bins, using the 20-60% centrality selection (see 
also 12511 ). Such a binning in pj and centrality was adopted in order to have a more significant 
comparison with previous preliminary results from the STAR Collaboration 12^1 . While STAR 
results are compatible with zero everywhere, a hint for a non-zero V2 can be appreciated at inter- 
mediate pj in the ALICE results. As pointed out above, partonic transport models indeed predict 
such a kind of behaviour. In Fig. |5jright), the forward rapidity inclusive J/0- V2 is shown as a 
function of centrality, in the range 1.5 < pj < 10 GeV/c. Again, clear hints for non-zero V2 can 
be appreciated for semi-central collisions. Depending on the adopted binning for the study, the 
significance for non-zero v'2 reaches about 3cr. 

Finally, the comparative study of the yields of various charmonium resonances in heavy-ion 
collisions is known to be an interesting tool to test theoretical models which implement differ- 
ent scenarios (sequential suppression, statistical production). Up to very recently the only results 
came from SPS energy, where a suppression of the iff(2S ) with respect to the 3/iff was observed by 
NA50 [27], increasing with centrality and with hints of a saturation for central events. Very re- 
cently, CMS studied the double ratio 1^(25 )/J/t^ between Pb-Pb andpp collisions B28I1 . observing 
values higher than 1 in the region 1.6 < \y\ < 2.4, 3 < pj < 30 GeV/c, although with rather large 
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Figure 5: Left: the pi dependence of the inclusive J/t/r V2 in 2.5 < y < 4, for Pb-Pb collisions (20-60% centrality). The 
results are compared to STAR values obtained at central rapidity, in the same centrality range. Calculations from partonic 
transport models are also shown. Right: the centrality dependence of the inclusive ]/if/ V2 in 2.5 < y < 4, 1.5 < pr < 10 
GeV/c. The bottom panel shows the average transverse momentum corresponding to each centrality bin. For both plots 
(left and right panels) the boxes around the points represent the systematic uncertainties. 



normalization errors, mainly due to the low integrated pp luminosity at -\fs — 2.76 TeV. ALICE 
can explore a contiguous range in rapidity (2.5 < y < 4) and also extend the pj reach of this 
measurement down to zero. In Fig.[6](left) the invariant mass spectrum for the region < pj < 3 
GeV/c, corresponding to a 40-60% selection in centrality, is presented. The result of a fit which 
includes a if/(2S ) component is also shown, and a clear hint for a signal can be seen, in spite of 
the rather unfavourable signal over background ratio. In Fig. [6] (right) we show the values of the 
double ratio if/(2S)/] /if/ measured by ALICE in two pi ranges. The pp values entering this ratio 
were calculated from the yfs = 7 TeV data sample, because of the too low integrated luminosity 
at -\fs = 2.76 TeV. The possible Vi'-dependence of the ratio was evaluated from a comparison 
of CDF, CMS and LHCb results and included (~ 15%) in the systematic uncertainty. Several 
other sources of systematic uncertainty cancel in the ratio, the main ones surviving being those 
on signal extraction and on the choice of the MC inputs for acceptance calculation. For central 
events, only an upper limit can be quoted for 3 < pj < 8 GeV/c, while background levels are too 
large to extract a meaningful result in the < pj < 3 GeV/c bin. Due to the size of the errors, 
no firm conclusion on an enhancement or suppression of tf/(2S ) / i / if/ in Pb-Pb with respect to pp 
can be drawn. However, our data tend to exclude a large enhancement in central collisions. 



5. Conclusions 

The ALICE Collaboration has performed a study of inclusive ifif/ production in pp and Pb-Pb 
collisions at the LHC, down to zero pj. The main results obtained include a suppression which 
saturates for semi-central to central events to values larger than the ones observed at RHIC. 
Differential studies of Raa vs centrality for various pj bins seem to favour a scenario where 
(re)combination processes play a sizeable role. The observed hints for a non-zero V2 signal 
at intermediate pi are in agreement with such a picture. Forthcoming p-A data will help to 
quantify cold nuclear matter effects and hopefully to sharpen the physics interpretation of the 
Pb-Pb results. 
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Figure 6: Left: the fit to the dimuon invariant mass spectrum for Pb-Pb collisions, for < pj < 3 GeV/c, in the centrality 
range 40-60%. Right: the centrality dependence of the double ratio [if/(2S)/]/if/]pb-piJ[if/(2S )/]/if/]p P of the i//(2S) and 
]/if/ yields at forward rapidity, compared to results from CMS (1.6 < \y\ < 2.4). The dotted lines refer to the pp reference 
uncertainties, while the shaded areas around the points represent the other systematic uncertainties. 
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